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Abstract 

A relativistic generalisation of a well-known method for approximating the dy- 
namics of topological defects in condensed matter is constructed, and applied to 
the evolution of domain walls in a cosmological context. It is shown that there 
are self-similar "scaling" solutions, for which one can in principle calculate many 
quantities of interest without recourse to numerical simulations. Here, the area 
density in the scaling regime is calculated in various backgrounds. Remarkably 
good agreement with numerical simulations is obtained. 
PACS numbers: 98.80.Cq, 11.27.-hd, 64.60.Cn 







Topological defects formed at a cosmological phase transition are one of the two 
known ways of creating scale-free primordial density perturbations [H, @] . However, unlike 
rival inflation-based models 0, analytic calculations have not yet made much impact. 
There do exist analytic treatments of the defects and Goldstone modes arising from the 
spontaneous breaking of global symmetries 0, which are exact in the limit of large 
N, where N is the the number of scalar fields, but they remain relatively undeveloped. 

In this letter a new technique is outlined which promises to form the basis of an exact 
dynamical theory for all topological defects described by a Nambu-Goto action. It is 
based on a method well-known in the condensed matter literature: the w-theory of Ohta, 
Jasnow and Kawasaki (OJK) and its descendents 0. In this approach, the defects 
are replaced by a multicomponent scalar field (the "u" field) specially designed so that its 
zeros track the positions of the defects. Although the equations of motion are non-linear, 
when combined with a Gaussian ansatz for the field probability distribution, they are 
susceptible to a mean field theory treatment. One can then calculate analytically many 
important quantities, such as the defect density and correlation functions, purely from 
the two-point correlation function of the w-field, which is itself calculable. 

To apply this theory to cosmic defects, all that is required is to make it relativistically 
covariant. The fictitious field that replaces the defects still has a Gaussian ansatz for its 
probability distribution function, and a self-consistent and self-similar solution for the 
linearised equations of motion can be found. With this in hand, one can calculate the 
defect density, using a generalisation of well-known techniques, although the presence 
of time derivatives of the w-field correlators complicates the procedure somewhat. In 
principle one could then go on to calculate more or less anything of interest: for example, 
in the cosmological context we would like to know the two-time correlation functions 
of the cosmic string stress-energy tensor p|. However, in this letter we shall content 
ourselves with calculating the scaling value for the area density of domain walls, in an 
approximation where one neglects the time derivatives of the w-field correlation function. 
The results for walls are encouraging when compared to the numerical simulations P, |10 



TT| . The results for strings are more difficult to obtain, and will be presented elsewhere 
Another feature of the theory is that it describes the behaviour of defects formed 



from initial conditions with a slight bias favouring one vacuum over another [|T2, |TT 
It is found that the defects disappear at a conformal time r]c = ''7i(f^^/(M^)c)^^^ where 
U = (u) is the initial bias in the field, and (u^)c the initial fluctuations around that value. 
Indeed, part of the motivation for this work was to account for this behaviour observed 



in some interesting simulations recently carried out by Larsson, Sarkar and White ||TT 



The essence of the technique is to replace the walls with a real scalar field, which 
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vanishes precisely at the coordinates of the walls X'^(o""), where fi takes the values 
0, . . . ,D and a the values 0, . . . , D — 1. Thus we begin with the equation u{X^) = 0. 
Differentiating once with respect to the world- volume coordinates cr", we find 

d^X^^OMX) = 0. (1) 

Thus the vector d^u{X) is a spacelike normal to the wall. 

The embedding of the p-brane in the background space-time induces a metric in 
its world volume, 7^/3 = g^u{.X)daX^^dpX'' , where g^^ is the space-time metric. Using 
the embedding metric we can covariantly differentiate (HD by acting with the operator 
(-7)-i/2^a(-7)i/2^"/3, to obtain 

nx'^d^u + -i'^^d^x^^dpX'd^d.u = o, (2) 

where □ is the covariant d'Alembertian and 7 = detja/B- The defect equations of motion 
are ii 

□x^ + r;;^7"^9„x^9^x" = 0, (3) 

where F^^ is the affine connection. In equation we can identify the tensor Pjl^ = 
'-y^'^do.X^dpXy, which is the tangential projector onto the wall. We can replace this by 
gij-v _ where = d^uduu/{duy. Using the defect equations of motion (||) we may 
then write 

This is the fundamental equation of motion for the field u which replaces the defects. It 
is actually unnecessary that this field has anything to do with the underlying Higgs fields, 
and so it may be called a "fictitious" field. In other approaches, which apply to global 
defects, the Higgs field is related to the fictitious field by the non-linear transformation 
{/)(a3,t) = /(u), where / solves the static defect field equations as a function of the 
transverse coordinate u [0. 

The equations of motion are not easy to solve, as they are non-linear. However, 
they have the distinct advantage over the original equations of motion (^) in that they 
are local: the defects may self-intersect and reconnect, a process which is highly non-local 
in its world volume. The approach taken in the condensed matter context is essentially 
that of mean field theory: bilinears in the field u are replaced by their averages. Thus 
one is assuming that the field is a Gaussian random field, and remains so throughout it 
evolution. This seems a good starting point for the relativistic version as well, although 
it should be borne in mind that the approximation is not well-controlled. 
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We begin the mean field theory manipulations by defining the basic equal-time two- 
point correlation function: 

(^(a;, r])u{x', r])) = C{\x- x'\,r)), (5) 

where the angle brackets denote an average over a spatially isotropic Gaussian probability 
distribution function. We shall also define M^jy to be the equal-time two-point correlator 
of duU: 

{d^u{x,r))d^u{x',r))) ^ M^^{\x - x'\,r)). (6) 
A two-point correlator with three derivatives will also be useful: 

{d^u{x,ri)d^dpu{x',ri)) = 'jf^^p^x - x'\,r)). (7) 

When referred to with no exphcit spatial variables, the correlators are to be taken at 
zero separation. In this coincident limit, the assumed spatial isotropy of the distribution 
function dictates their forms. The non-zero components of M^^^, are 

Moo = T{i]), M^, = S{r])6^n, (8) 
where m,n = l,...,D, and the non-zero components of 7^^,^^ are 

7000(^) = lOmniv) = -^S{r])6mn, lmOn{v) = ^S{r])Smn- (9) 

We now linearise the equations of motion by taking the Gaussian average, and then find 
a self-consistent solution for the fields u{x,r]). We need the following identities: 

{{dufd^d^u) = Mdpd,u + 2^pp,gP''d^u (10) 

{{du^dpu) = MdpU + 2gP''Mr^pd^u, (11) 

where M = Mp^^g^"". 

In a flat Priedmann-Robertson- Walker space-time, the affine connection is — 
9iJ.i^9^){^/'^)- We can now see that the hnearised equations of motion have 

the form 

ii+^u-v^W^^O, (12) 

T) 

where n{r]) and v depend on D. For Priedmann models, one can show that 

l,(r)) = -2n{S/S) + a{v)D[l-3{T/S)], (13) 
= [D-l-iT/S)]/D, (14) 

where a{ri) = rja/a. 
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In a scaling solution we would expect S and T to have a power-law behaviour with rj. 
Thus, so as long as we are not near a transition in the equation of state of the Universe, 
yU and are constant. Imposing the boundary condition that u be regular as — » 0, 



(O) then has the simple solution 



where Ai, = 2'^r(z/ + l)uj^{rii), and (1 — fiy/4: = . The form of the initial power 
spectrum Pi(A;) = \u\Jk^\)^ is taken to be white noise, which ensures that the field is 
spatially uncorrelated to begin with. 

We may now evaluate T jS and f^, and implicitly solve for v. Firstly, we must decide 
the sign of v. It turns out that it is inconsistent to take (1 — /i)/2 = v (one way of 
showing this is to compare C{j]) calculated by explicit differentiation, and by evaluation 
of {uii)). Thus C scales as rj^^ ., S and T as r/~(^+2). Using standard integrals of Bessel 
functions, and defining the parameter /? = 2i/ — D — 1, we find 



T _ {D + 2){D-2) 
S ~ 3{D + 2) + 2/3 



(16) 



provided /? > 0, so that the integrals for 5* and T are defined. To find P, we solve the 
equation 

/i = /5 + D + 2 = aD[l- 3{T/S)] + 2{D + 2). (17) 

The radiation-dominated {a = 1) and matter-dominated {a = 2) Friedmann models 
require a certain amount of algebra: however, Minkowski space {a = 0) has the simple 
solution p = D + 2, giving T/S = (D - 2)/5 and = {W - 3)/5D. The smallness 
of T/S will help us in the next part of the argument, where neglecting T (and another 
correlator involving time derivatives) in comparison to S will make calculations much 
easier. 

Armed with the "mean-field" solution for u{x, rj) we can now calculate anything that 
can be expressed in terms of local functions of the field and its derivatives, provided of 
course that we are able to perform the Gaussian integrals involved. In this letter we 
content ourselves merely with evaluating the area density, which is given by 



Making the coordinate transformation from to (a", u) near the wall, this can be 
rewritten as 

A = 6{u)\du\. (19) 
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Table 1 : The calculated values of the comoving domain wall area density in the self-consistent 
linearised solution for u. Values listed are for Minkowksi space (a = 0), and radiation and 
matter-dominated FRW models (a = 1,2 respectively). 



In order to calculate the Gaussian average of A, one takes Fourier transforms 

V{D/2) f dk f d^^'q I d d 



^■^^ 27ri+^/2 j 2txJ (g2)i3/2 \dq' dq ""^^ / ' 

One potential difficulty is that q^ can vanish, due to the Lorentzian metric. The equation 
must therefore formally be defined by analytic continuation to imaginary time, and an 
accompanying Wick rotation of the time components of the Fourier transform variable 
q^'. 

We neglect terms of order T/S, which complicate the calculation considerably. With 
less justification, we shall also neglect terms involving jCS which appear in the full 
calculation, and are of order 0.5, as one can check from the formula (^) for SjC below. 
The averaged comoving defect area density A{j]) is then just 



(22) 



This is identical to condensed matter results for domain walls [TB|, T\ 
One can show that 

S_ ^ 1 (/? + l)(/3+l + D/2) 
C~ if 2Pv^ 

and we immediately see that the correct scaling behaviour for the wall area density 
{A{r]) oc is reproduced. Furthermore, we are able to compute the coefficient, and 
compare with numerical simulations. The comoving area densities for walls in two- and 
three-dimensional Minkowski space, radiation- and matter-dominated Friedmann models 
are displayed in Table 0. 

The theory seems to work well: scalar field theory simulations in D = 3 give a 
comoving area density of approximately 1.5/ rj in both radiation |T0| , pd]| and matter [0 
eras, which sits nicely with the calculated values. However, the agreement is partly 
fortuitous, as there are probably significant errors in both the theoretical and numerical 
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values. In D = 2 the agreement is also good, although there is evidence for a small 
deviation from the r]~^ behaviour in = 2 P, [TO . 



One may also ask how the network behaves when a small bias is introduced into the 
initial conditions, that is, if {u{x,rii)) = U. In numerical simulations of domain walls 



pHI , |TT|, it is found that even for very small initial biases, for which the walls percolate. 



the system still evolves away from the percolating state and eventually the large walls 
break up and disappear. Similar behaviour is well-known in in the study of quenches of 
condensed matter systems with a non-conserved order parameter P, |13|, |15|, p!6| . 

The theoretical description of this behaviour is fairly straightforward. Introducing 
a bias into the initial conditions alters the Gaussian average of ( pOl) to {A{ri))i/ = 
{A{r])) exp {—U^ / 2C {t])) , where {A{ri)) is the zero bias result from (|2T|). If the sys- 
tem is close to being self-similar at some initial time rji when the magnitude of the bias is 
U and the fluctuation around that value is C{rii), one can calculate the time t]c at which 
the defect density falls to a fraction of its scaling value to be 

V, = Vi{uV2C{ndy'^''. (23) 

Recent simulations by Larsson, Sarkar and White are consistent with the above calcula- 
tions of {A{ri))u and ?7c in D = 2, but do not have sufficiently good statistics in Z) = 3 
to be able to check the results |]TT|. Coulson et al. [T^ did not attempt fits of the correct 
form, A oc ?7~^ exp(— (t^/t^c)'^), to their simulations, although they did note that the walls 
disappeared faster than a simple exponential in D = 3. 

To summarise, this paper outlines a new analytic technique for describing the dy- 
namics of topological defects after a cosmological phase transition. It is a relativistic 
version of a well-known approach in condensed matter physics, due to Ohta, Jasnow 
and Kawasaki 0. The scaling area density of domain walls in three spatial dimensions 
is calculated, and agrees quantitatively with numerical simulations of a real scalar field 
@, 10, |ll|]. Further work is clearly needed, particularly on strings, but it already seems 
clear that developments of this technique will be of great benefit to the study of cosmic 
defects. 

I am deeply indebted to Sebastian Larsson, Subir Sarkar, and Peter White for com- 
municating their results to me in advance of publication, and for useful discussions. 
The seeds of this work lie in the Topological Defects programme at the Isaac New- 
ton Institute, Cambridge, July-Dec 1994. I am also grateful for the hospitality of the 
NORDITA/Uppsala Astro-Particle Workshop 1996, where this work was completed. The 
author is supported by PPARC Advanced Fellowship number B/93/AF/1642, by PPARC 
grant GR/K55967, and by the European Commission under the Human Capital and Mo- 
bility programme, contract no. CHRX-CT94-0423. 
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